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Abstract 



A new measure of the crystal-field strength, complementary to the conventional one, is defined. 
It is based on the rotational invariants |-Bfco| av or Bko\ av , k = 2,4, 6, of the crystal- (ligand)-field 
(CF) Hamiltonian 7Ycf parametrizations, i.e. on the axial CF parameters modules averaged over all 

reference frame orientations. They turn out to be equal to TLqI and |HcF| av , respectively. While 

av 

the traditional measure is established on the parametrization modules or on the second moment of 
the CF energy levels, the introduced scale employs rather the first moment of the energy modules 
and has better resolving power. The new scale is able to differentiate the strength of various 
iso-modular parametrizations according to the classes of rotationally equivalent parametrizations. 
Using both the compatible CF strength measures one may draw more accurate conclusions about 
^ . the Stark levels arrays and particularly their total splitting magnitudes. 

G\ : PACS: 71.70.Ch 
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1. Introduction 



Solid state experimentalists, especially spectroscopists, still need a reliable scale quantitatively char- 
acterizing the effect of crystal-field interaction, i.e. defining the so-called crystal-field strength. Such a 
parameter could directly verify and compare various parametrizations of the crystal-field Hamiltonian 
TCcf, which may come from different fittings experimental data when the orientations of reference 
frames associated with these parametrizations are unknown in the majority of cases. 

Although such a conventional scale for measuring the strength of the crystal-field has been already 
introduced over twenty years ago [1,2], in some cases it seems to be insufficiently precise. It employs the 
basic rotational invariants of the TCcf, i-e. the modules of its 2 -pole components TC C F, defined as Mk = 

[J2q \Bkq\ ) , as well as uses the global TCcf modulus M = iJ2k J2 q \^kq\ ) ■ in the first case the 

/ , \l/2 

partial crystal-field strength is defined as Sk = ( ) Mk, while in the second case the global crystal- 
field strength is given by S = (Efc^) 1/2 . Throughout the paper the tensor (Wybourne) notation 
for the crystal-field Hamiltonian and the crystal-field parameters (CFPs), TCcf = J2kJ2 q BkqCq k \ is 
consistently used [3]. The summations over k and q indices run, in each individual case, over strictly 
specified values according to the kind of central ion and its point symmetry. 

Both the parameters Sk and S themselves are not a direct measure of the real magnitude of 
the initial state splitting, since the crystal-field effect depends also on the properties of an object 
(a paramagnetic ion) upon which the TCcf acts. Namely, the response of the system to the TCcf 
perturbation reflects the symmetry of the electron density distribution of the central ion open-shell. 
For instance, an S"-type ion like Gd 3+ feels no crystal field (in the first order of perturbation) no 
matter how strong is the surrounding field. 

The effect of splitting can be most simply expressed by the so-called second moments a\ or a 2 of 
the CF sublevels within the initial state upon switching on the TCq F (or TCcf) perturbation. In fact, 
the second moment is easily represented by the scalar crystal-field strength parameters, either Sk or S 
(section 2). However, although the effective TCq F multipoles (for k = 2,4,6) contribute to the energy 
of individual Stark levels independently (as an algebraic sum), the simple linear relations between a\ 
(or a 2 ), and S 2 (or S 2 ) are always fulfilled. As it is proved these relations strongly confine both the 
maximal (A£ max ) and minimal (A£ m j n ) nominally allowed splittings of the initial state (section 3). 
Moreover the actual crystal field splittings AE can be additionally restricted (section 5). Naturally, 
all the iso-modular TCcf parametrizations correspond to the same crystal-field strengths Sk and S. 

(k) 

However, apart from the modules Mk and M, there exist also other rotational invariants of the TC CF 
or TCcf which distinguish the whole classes of the rotationally equivalent TCcf parameterizations, in 
other words the parameterizations referring to the same real crystal-field potential, but expressed in 
variously oriented reference frame. Interestingly, the new invariants turn out to be the average values 
of the axial parameter modulus |-Bjto|av, k = 2,4,6, in the case of TCq F , or | Efc-^fcolav for the global 
TCcf obtained after the averaging over all orientations of the reference frame, i.e. over the solid angle 
4ir. As it is shown in the paper the average value of the axial parameter modulus or the average of 
the modulus of their sum are just equal to |H^p| a v and | Hcf lav, respectively (section 4). 

The new scale of the crystal-field strength based on the above invariants is in principle consistent 
with the conventional one but it reveals more resolving power. Applying the new measure to the 
iso-modular parametrizations may lead to different strength parameters what is exemplified below 
for several cases (section 5). The introduced more subtle strength gradation established rather on 
the first moment of the sublevel energy modules gives, comparing to the second moment, additional 
information about the Stark levels array for various iso-modular TCcfs, including the magnitude of 
the total splitting gap of the states. In this paper we confine ourselves to the pure model states 
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of the zero-order approximation with a well defined angular momentum quantum number and the 
corresponding degeneration. These could be for instance Russell- Saunders coupled states \aLSJ) 
coming from the 2S+l L terms, where a stands for the remaining quantum numbers needed for their 
complete determination. Such states have a well defined quantum number J and the degeneration 
2J + 1. The derivation of the analogical expressions including J-mixing effects [4] or a transformation 
to other functional bases of the zero-order approximation can be accomplished by using standard 
angular momentum re-coupling techniques [4-8]. In section 5 we analyse by way of example the 
crystal-field splitting of p 1 , d 1 and f 1 one-electron configurations and a typical complex state 3 #4 

(k) 

for various iso-modular H C p, k = 2,4,6. In the first three cases we avoid complex states re-coupling 
procedure which is a side issue to the problem under consideration. Since we study the differentiation 

(k) 

of the effects due to various iso-modular Hamiltonians TC C p, all CFPs values along with the Stark 
levels energies are given in units. 



2. Conventional definition of the crystal-field strength parameter 

The comparison and scaling of the crystal-field impact can be based upon the two types of scalar 
quantities, and/or M, since both of them are rotationally invariant. A scalar crystal-field strength 
parameter of this kind was given firstly by Auzel and Malta [1,2] as (in original notation): 



N v = 



k,q 



2ir 



2k + 1 



1/2 



which is nothing more but M in the space spanned by spherical harmonics Y q k . In other words, N v 
is a norm representing a distance in the space. Currently there are two definitions widely used in the 
literature [9-12]: 



1/2 



Sk — 



T x>. 



1 



2k + 1 Y ' J \2k + l 

for k = 2, 4 and 6 in the case of 2 fc -pole Hcf component and 

,V2 



B 2 k0 + 2j2(ReB kq ) 2 + (ImB kq ) 2 

q>0 



1/2 



(1) 



s = (si + sj + s 2 ) 



or 



S 



U s 2 + s 2 + s 2 ) 



1/2 



(2) 



for the global Hcf = Z)fc^cF Kl- ^ word of caution seems to be worthy at this point. Namely, the 
values of Sk or S can differ according to the type of the Hcf parametrization (operators) applied. 
They can be compared with each other only after proper recalculation. Since both these quantities are 
independent of the assumed axis system they allow to check whether the original CFP data sets and 
the transformed ones are compatible. The strengths Sk or S enable also a broad comparison of CFP 
data sets when the axis systems have not been explicitly defined, and undoubtedly they play a central 
role in the CF theory. What is also important and useful they are linked to the second moment of the 
Stark levels within a particular initial state \aSLJ) [4,13]. 

The second moment of the sub levels \n) within \aSLJ) state upon introduction of a TCcf pertur- 
bation is defined by 



a 2 (\aSLJ)) 



2J + 



-Y,[En-E{\*SLJ))\ 



(3) 
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Hence 



where the center of gravity of the Stark levels belonging to the state \aSLJ) is given by 

E(\<*SLJ)) = ^J-_££ n , 

and E n is the |n) sublevel energy. Since TLcy is diagonal in the |n) basis and the second order effect 
of Tic? interaction is neglected [4] 

a 2 (\aSLJ)) = ^L_Tr{H2 F }. 

a 2 (\aSLJ)) = -^^sl((aSLJ\\CM\\aSLJ)) 2 , (4) 

what implies from the orthogonality of 3-j symbols [2,5,13]. The symbols {aSLJ\ \C^ \ \aSLJ) are 
the double-bar or reduced matrix elements of the spherical tensor operators. According to Wigner- 
Eckart theorem [14] they are independent of the reference frame orientation. Their origin and physical 
meaning stem from the following relationships [5-8,15]: 

(aSLJMj\C^\aSL'j'M'j) = (-l) J " Mj ( _^ J ^ ) (aSLJMj\\C^\\aSL'J'}, 

where the reduced matrix element follows the 3-j factor. Further use of tensor formalism yields 

(aSLJ\\cW\\aSL'J') = (-l) s+L '+ J+fe [(2J + 1)(2J' + 1)] 1/2 | ^ ^ * j (aSL\\C^\\aSL'), 

where the double-reduced matrix element follows now the 6-j symbol. We can also pass to the matrix 
elements of the unit operators [5-8,15], i.e. normalized equivalents of C^ k \ since 

(aSL\\C (k) \\aSL'} = (aSL\\U (k) \\aSL'}(l\\C (k) \\l), 

and / is the angular momentum quantum number of the open-shell electrons. The reduced matrix 
elements of the operators have been compiled by Nielson and Koster [16], whereas the 3-j and 
6-j symbols can be found in the tables by Rotenberg et al [7]. 

The simple relation between the a\ and S\ (Eq.4) can be also proved employing Vieta's formulas 
for roots of the Hq^ matrix characteristic polynomial 

E n + ai E n ^ + a 2 E n - 2 + . . . a n = 0, 

which is here of order of n = 2 J + 1. All its coefficients and roots must be real what follows obvi- 
ously from the TtcF hermiticity. Interestingly, some characteristics of the sublevels spectrum may be 
described in terms of the elementary algebra. Firstly, as the energy center of gravity of the initial 
state is conserved, i.e. (Ya=i E% = 0) 2 > the a\ coefficient standing at E 11 " 1 must vanish. Next, since 
= {Ya Ei = 0) 2 = Ya=i E 2 + 2J2i>j EiEj, the second moment, i.e. the sum of the root squares 
(divided by 2 J + 1) is equal to 2j+\ ■ ^ can ^ e a ^ so sri ° wn that 

" 2 ° 2 = ^T M ^ J ll c(fc) H J ) 2 ' 
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where the simplified notation for the reduced matrix element representing only the last quantum 
number has been introduced. Hence, between a\ and S 2 a simple formula holds (Eq.4) 

4 = 2jVi^< j h cW h j > 2 - (5) 

In other words, is proportional to S^. Finally, a free term of the characteristic polynomial is given 

as a n = (— l) n EiE2 . ■ ■ E n , what may be helpful analyzing the solutions. For instance, if one root 

equals zero then a free term vanishes. 

The problem becomes more complex for the global crystal-field strength S (Eq.2), since then the 

components S% contribute to the sum with their weights (J\\C^\\J) 2 (Eq.4). This is why there is 

no straightforward proportionality between a 2 and S 2 in this case. Nevertheless a is a positively 

defined quadratic form of S\. and, in consequence, the inputs of particular Hcf 2 fc -poles into a 2 

cannot compensate themselves. The condition that a 2 is constant for various iso-modular Hqf does 

not exclude, however, a possible differentiation of the CF sublevels sequence and structure, as well 

as the initial state total splitting. In fact, a and S could be correlated similarly as and in the 

previous case, only if the elements (J||C^ fc ^|| J) were equal for all k. Nevertheless, the second moment 

of the Stark levels within a particular state \aSLJ) is simply given in terms of S2, S4 and Auzel 

and Malta [2] made an attempt to average the a 2 quadratic form by bringing down its respective 

r 1 1/3 

ellipsoid J2k S%{J\\CW\\J) 2 in the fc-space to a sphere of the same volume J2k S% U k {J\\C^\\J) 2 
and having a radius equal to the geometric mean of the three ellipsoid axes. In practice, unfortunately, 
this elegant approach does not always lead to acceptable results. In the literature the overall effect 
of the crystal-field interaction is often characterized by a quantitative comparison of the crystal-field 
strength [17-20]. Additionally, a systematic correlation between the free ion parameters and the CF 
strength is observed, namely increase of the crystal-field interaction results in the reduction of the 
free- ion parameters [17]. The CF strength increases in the RE series with decreasing ionic radius of 
the RE 3+ host cation [19]. The physical meaning of the CF strength scalar parameter is also supported 
by the fact that it rises with pressure applied to a sample [21,22]. The CF strength parameter has 
also been used to compare the root mean square error obtained for crystal fields of different strength. 
However, its use in such a case is restricted only to comparisons of the identical site symmetries 
[10,23]. Furthermore, within the approximation to the second order in the crystal-field, the shift in 
the center of gravity of a particular 2S+1 Lj state due to J-mixing effects is a simple linear function of 
the Sf. [4,17]. The concept of the Sf, or S can be extended to define the quantities Ck and Cq [13] as 
normalized "scalar products" of any two compared parametrizations. These quantities represent the 
"angles" between the two considered parametrizations and are a convenient measure of the closeness, 
i.e. the correlation of any two CFPs sets. 

3. The correspondence of the Stark levels second moment of | J) state to its nomi- 
nally allowed splittings 

The second moment of CF levels, a 2 , essentially limits a formally allowed range of the initial state | J) 
total splittings AS for different but iso-modular Hqfs- Such energy splitting confinement differs for 
non-Kramers and Kramers ions what is specified in details below. 

Let us firstly study the case of any integer J, i.e. non-Kramers ions. Having to keep a constant a 2 
the minimal hypothetical splitting A£ m ; n of the (2 J+l)-fold degenerate state takes place when J levels 
assume identical energy of jj^A£ m i n , and the remaining J + 1 levels take the energy A£mim or 
vice versa. Further this is referred as Type I splitting. Then 

2 J(J + l)(A£ min ) 2 2J + 1 
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In turn, the maximal hypothetical splitting A£ max occurs for one level of A£ max /2 energy, one of 
— A£ max /2, and the rest (2J — 1) levels with zero energy. Further this is referred as Type II splitting. 



Then 
and hence 



2 = ^IjTi^ 2 ' A£ max = ^2(2J + l), (7) 



A£ max / 2J(J + 1) 

A£ min V 2 J + 1 ' 
Let us also consider, following Auzel and Malta [2], the case of the homogenous splitting A£hom> when 



2 2(1 + 4 + . . . + J 2 ) (A£ hom /2J) 2 / 3J 

Below it will be referred as Type III splitting. For example, if J = 4 then AS m - m = a ~^j^ = 2.01c, 
A£ max = a3^2a = 4.24a, A£ hom = AaJ3/5 = 3.10a, and finally the ratio 45™ = 2.11. 

For J = 1, A£ min = a^, A£ max = A£ hom = ay/E. The ratio fj^ = 2/^3 = 1.16 and this 
narrow interval strictly limits the AS variation. It has a simple graphical interpretation. As is known, 
three real roots of a third order equation must fulfill the conditions (Cardan's formulas) presented 
in Fig. la, where the angle (p is a function of the equation coefficients. The maximal and minimal 
splittings AS correspond to the solutions shown in Figs lb and lc, respectively. 

Let us now pass to the Kramers ions with a half integer J. Here, two cases should be analyzed. 
Firstly, if an even number of doublets (2 J + l)/2 occurs the minimal | J) state splitting takes place 
when (2J + l)/4 doublets have the energy AS m - m /2, and the next (2J + l)/4 doublets the energy 
-A£ mm /2. Then, 

2 4[(2J + l)/4](Ag min /2) 2 ac 0fT 

In turn, the maximal splitting, A£ max , will appear if one of the doublets will be of energy A£ max /2, 
and the second one of energy — A£ max /2 with all the rest of levels with zero energy. This time 

2 4(A£ max /2) 2 ac „ /oTXT 

= 2 J + 1 ' A£ max = oV2J + 1, 

and therefore now 



Af max _ V2JTT 
AS ■ 2 

Secondly, for an odd number of Kramers doublets A£ m j n = 2^(2 J+i) an( ^ 

y (2J+3)(2J— 1) 



max 



aV2JTT, with = \^m^- 

The homogenous splitting AS\ iQm for an even number of doublets (J = (4k + 3)/2) and for an 

odd number of doublets (J = (4k + l)/2), where k = 0, 1, . . ., amounts correspondingly to cr^ 3 *- 2 



and 2cry 3 ^ 2 y +3 1 - > . By way of example, if J = 9/2 (i.e. for five doublets), then A£ mm = cr-j^ = 2.04cr, 

A£ max = a\/l0 = 3.16<7, and A£h om = 2.83cr. In the case of Kramers ions the AS variation range 
turns out to be smaller than that for non-Kramers ions, which is seen comparing the A£ m i n and A£ max 
for J = 4 and J = 9/2. Finally, taking the most extreme case of J = 15/2 for /-electron configurations 
(e.g. for Dy 3+ , Er 3+ ) with eight doublets, we would obtain AS m i n = = 2.02a, A£ max = 4a, and 
AS hom = 2.37a. 
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4. The new scale of the crystal-field strength. Comparison of both 
the scales S k = ^M k and S' k = 



av 



4.1. Average values of the axial parameter modules |-Bfc | av an d | J2k -^fcoUv? where 
k = 2,4,6 — the rotational invariants of the equivalent T~Ccf parametrizations 

Rotating the reference frame by the two Euler angles a and (5 we obtain all the equivalent Wcf 
parametrizations (with the accuracy to the third Euler angle 7 about the z axis) [5,24]. Their axial 
parameters for a 2 fc -pole component are given as: 



Bi 



k k 

q=—k q=—k 
k 

= C^\f3)B k0 + 2 £ Cf\f3) cos q(a + <p q )\B kq \, 



where T>Q^(a, /?, 7) are the middle row rotation matrix elements, Cq k \f3,a) 



the spherical tensors, whereas Cq k \(3) = (— l) 9 ^ P^(cos /?), and P k (cos (5) are the associated 

Legendre functions, -B^ = l-B^le 19 ^ 9 , and -Bfc_q = (— l^l-B^e -19 ^ 9 . The primed parameters corre- 
spond to the transformed parametrization while the unprimed to the initial one. It can be directly 
proved that (B' k0 ) av = and (J2 k ^fco) av = ^> wnne the average absolute values 

2n 7T 



I 4tt 
^2fc+l 



1/2 



(10) 



Yq(j3, a) are 



\ B 'koL = -tj f\B' k0 \(a,P)\ sin Pd(3da 



AO 





2-7T 7T 



4tt// 




sin/3d/?<ia , 



(11) 



as the rotational group invariants are discriminants of the equivalent parametrizations classes [24] . By 
the mean value we understand the magnitude averaged over all possible orientations of the reference 
frame, i.e. over the solid angle Air. Interestingly, they can be used to estimate the CF strength 
independently of the parametrization modulus. 



4.2. Average values of the modules 

Since the expression 



n-f(k) 

n CF 



and \Hcf\ 



H 



(fc) 

CF 



£ B kq cW(P,a), 



where (3 and a are the spherical angle coordinates in the central-ion reference system, is identical with 
that for B' k0 (Eq.10), the following important identity holds 



2-7T 7T 



Ti 



(k) 

CF 



^/ / \n§\(a,(3)\ am fidfida = \B' k0 \ 




(12) 
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The average value of the modulus of the 2 fe -pole Hq^ component turns out to be equal to the average 
value of the modulus of the relevant axial parameter B^q. This identity, Eq.12, obvious when we prop- 
erly interpret the rotation angles in both cases of averaging, associates l-B^oLv wu -h the complementary 
measure of the CF strength S' k : 



S'k — 



n CF 



av 



= l^olav (13) 



Although the expression for S' k in the above form is limited to a given 2 fc -pole Ti-Q-p component, it may 
be generalized for the global Hcf 



S' 



CF 



k 



(14) 



Contrary to the conventional CF strengths Sk and S (Eqs 1,2), which are constant for all the iso- 
modular parametrizations, the new strengths S' k and S' calculated for the constant modulus (and 
modules) change their magnitudes within certain ranges discussed in the next section. To compare 
both the measures it is convenient to express S' k in the product form f k ■ M k , where fa is a specified 



factor. Now, these two measures will be compatible if the factor fa is close to y 2fc+1 , Eq.l, i.e. to 
0.447, 0.333 and 0.277 for k = 2, 4 and 6, respectively. This compatibility is demonstrated in the next 
section, where a thorough discussion of the relation between both the CF strength scales is provided, 
by way of example of the CF splitting of p 1 , d 1 and f 1 electron configurations with the spin-orbit 
coupling deliberately neglected , and the 3 #4 state for various iso-modular Ttg F s. 

5. Computational results and discussion 

5.1. Crystal-field splitting of p 1 , d 1 and f 1 electron configurations for various iso- 
modular Hq F s, k = 2,4, 6. 

We consider the model results of interaction of any iso-modular Hq F (k = 2,4,6) with M k = 1, on 
the initial states with well defined angular momentum quantum numbers. The magnitudes of all 
quantities under discussion, i.e. the new CF strength parameters S' k , the total splittings AE^ k \ the 

second moments a\ of CF levels and the averages of the absolute values of the Stark level energies 

(k) 

\En |av are given in units. Tables 1, 2 and 3 present a comprehensive review of S k values for 

(k) 

various iso-modular 7i CF s, with k = 2,4,6, respectively. Correspondingly, these five, ten and eleven 

(k) 

Ti-QpS compiled in Tables are the representative ones including those with the highest and lowest S k 

(k) 

values found during the survey. No other T~C CF s seem to yield S k out of these ranges. The strength 
parameters S' k change themselves within the rather narrow intervals: 0.368 — 0.385, 0.251 — 0.287 and 
0.195 — 0.239, while the relevant are constant and equal to 0.447, 0.333 and 0.277 for k = 2,4 and 

(k) 

6, respectively. The maximal S k parameters refer to the purely axial H-Qps when B k QS achieve 1. For 
other parametrizations this maximal value of 1 is not achieved in any reference frame. 

As implies from Tables 4, 5 and 6 there is a certain mapping between the above S' k ranges and 
the referring to them intervals of AE^> and \E^\ av . As is shown in the paper this quantitative 

(k) 

mapping is determined by the roots of the TC CF matrix characteristic polynomial, and the key part 

( J k J 1 \ 

of the matrix elements is the product {—l) Mj B kq , . The remaining factors coming 

y — iWj q Mj J 

into the matrix elements are common and play the role of a scaling factor. In the below examples 
concerning the CF splitting of one-electron states with J = J' = I for I = 1, 2 and 3, the role of such 
a scaling factor play the double-bar matrix elements (Z||C^||Z). 
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It should be pointed out, however, that the mappings S' k < — ► AE^ k \ S' k < — ► |-En^|avj AE^ <— 

(k) 

—> \En |av are neither straightforward nor explicit. With the increase of the initial state degeneration 
2 J + 1 they become less clear due to a big variety of possible splitting schemes. Nevertheless, one 
may presume a dominant tendency: the greater S' k the greater |.E^|av and the lesser AE^ (Type 

I splittings). In the reverse case, i.e. for a small S' k , Type II splittings are expected. However, such 

(k) 

reasoning does not take into account the unique characteristics of the Hamiltonian averages TC C t? , 

av 

and the space density distribution of unpaired electrons in the states of various J. From this point of 
view the analysis of Tables 4, 5 and 6 seems to be instructive, indeed. 

On the other hand, the allowed spans of the A£^ for a fix M k , i.e. <Tfc are known. In the light 
of the above mapping it turns out that not all of these values A£^ k \ and corresponding to them 
splitting schemes, can actually occur. Namely, depending on the initial state quantum number J and 
the multipole's rank k some specified limitations of the AE^ are observed (Tables 4, 5 and 6). They 
are listed briefly below. 

For I = 1 (p-electron) and k = 2 the full nominal range of the A<?( 2 ) and all splittings of Types I, 

II and III are admitted. More particularly, AE^> can vary from 0.600 M 2 to 0.693 M 2 (Table 4). 

For I = 2 (d-electron) and k = 2 the magnitude of AE^> is constant and equals 0.572 M 2 for each 

(2) 

iso-modular T~C C p what corresponds to splittings similar to those of Type I. Other splittings, including 

(2) 

e.g. A£^ ' are impossible in this case (Table 5). 

For I = 2 (p-electron) and k = 4 again the full nominal range of the A£^ is allowed beginning from 
the smallest 0.363M 4 for the cubic ftgj., up to the biggest 0.564M 4 for HqI = ^B u c[ 4) + ^B 4 „ 4 C { fl 
(Table 5). 

For I = 3 (/-electron) and the value of k = 2 AE^ weakly depends on S' 2 , varying in all its range 

(21 

merely from O.6OOM2 to O.6O8M2, i.e. somewhat below the A£^J m (Type III splittings) (Table 6). 

Next, for / = 3 and k = 4, the possible AE^ varies within the range from O.358M4 to 0.482M4, 
i.e. around the A£^J m (Table 6). 

Finally, for / = 3 and k = 6 the allowed AE^ varies from 0.326M 6 to 0.501M 6 covering the 
majority of the nominal range together with its upper limit, but excluding the smallest splittings 
(Table 6). 

The obtained results may be generalized for states with J or L equal to 1, 2 or 3, multiplying 
A£( fe ) and by the scaling factors (J\\C^\\J) or (L\\C^\\L). 

5.2. Crystal-field splitting of 3 H 4 state in various iso-modular Hq^s, k = 2,4, 6. 

Let us end up with the analysis of splitting of nine-fold degenerate 3 H 4 state subjected to the iso- 
modular HqpS enclosed in Tables 1, 2 and 3. Table 7 shows the correlation between S' k , AE^ and 
|^ fc) | av . The scaling factors (J = 4\\C^\\J = A) = (J = 4\\U^\\J = 4)(f\\C^\\f), required here due 
to the coupled initial state (L = 5, S = 1, J = 4), are equal to —1.2365, —0.7389 and 0.7706 for 
k = 2,4,6, respectively. Hence a 2 = 0.184M 2 , a 4 = 0.082M 4 and a 6 = 0.071M 6 , while the global 
second moment of the Stark levels takes the form 



-(-1.2365) 2 Mf + -(-0.7389) 2 Mf + — (0.7706) 2 M 6 2 
5 9 13 



The ranges of the formally allowed A£^ corresponding to the above second moments o 2 ,o 4 ,<J§ are 
marked in Fig.2 by the solid lines. 

We can see in Fig.2 that from the set of all potentially allowed total splittings A£^ only certain 
AE^ may be realized (those between the dashed lines), and consequently, only certain splitting 
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schemes (roughly between Types I and III) may occur. For instance, in the case of all the three 
effective multipoles neither ASmlx nor A£^ n are possible, while A^oL can appear solely in the case 
of 2 6 -pole. Based on Table 7 it is seen also that for all the three effective H^-pS the biggest AE^ are 
achieved for intermediate Si values. 



6. Conclusions 

The conventional scales Sk or S with the associated second moments of the CF levels, <r fc or a, do 

(k) 

not distinguish the iso-modular TC C p or TCcf parametrizations, which, however, can be differentiated 
by an another scale - the spherically averaged S' k = Hqp and S' = |HcF| av - It is proved that 

clV 

the S' k variation ranges for all the iso-modular parametrizations are limited and lie slightly below 
the relevant Sk magnitudes. The span of these ranges amounts to 5, 10 and 20% of their values 
for k = 2, 4 and 6, respectively. There exists a direct mapping of S' k ranges into the total splitting 
AE^ ranges and \En\ecv intervals, which may be interpreted more clearly for the initial states with 
low degeneration. Such mapping allows to estimate the total splittings AE^ or AE to be expected 
and characterize their spectrum. It is shown that not all the nominally admitted total A£^ or AS 
splittings determined by the modules or M, i.e. the second moments cr^ or a, can actually occur. 
This essentially confines the set of the allowed splitting schemes. 
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Table 1: The spherical averages of five representative iso-modular 7i C pS, S 2 = \Ti-Qp\avi acc - to Eqs 
11-13, expressed in M 2 units. Only B 2q CFPs are given, B 2 - q = {-^) q B 2q 



No. composition S 2 = |H^p| av 



B20 B 21 B 22 



1 1 0.385 

2 7S 7S "75 °- 381 

3 75 7S "75 °" 374 

4 75 7S eW4 71 °- 369 

5 -4 0.368 



V2_ 
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Table 2: The spherical averages of ten representative iso-modular TCq^s, S 4 = |H^p| a „, acc. to Eqs 
11-13, expressed in M4 units. Only B± q CFPs are given, B^ q = (— \) q Bl q 

No. ^cf composition S' 4 = |Wcp| av 
B40 B41 B42 -B43 B44 

1 1 0.287 

2 iyf 0.280 

3 III eW2 J I °- 277 

4 - — - - - - 276 
^ 3 3 3 3. 3 u./wu 

5 I I I ~l e1 ^ 2 I °- 273 

6 0^00 0.269 

7 ^ 0.266 



3 3^ 3 3 



0.265 



Q I I I I I 

3 3 3 3 3 

9 ± ±e i7r / 4 i i i 0.261 



10 -i 0.251 

i/2 
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Table 3: The spherical averages of eleven representative iso-modnlar W g> B| S' 6 = \H^\ av , ace. to 
Eqs 11-13, expressed in M 6 units. Only B 6q CFPs are given, B 6 - q = (—l) q Bl q 

No. H-cf composition S' e = |H^p| av 

-^60 Bqi Bq2 Bqz B e4: Bq 5 Bqq 

"1 1 0239 

2 -J_ -J- -J- .J- L_ _J_ ^_ 231 

Vl3 Vl3 VlZ VlZ VlZ VlZ 

3 -U -U — i- -i. 4- -r- -r- °- 228 

v 7 !! Viz Viz Viz Viz Viz Viz 

A _J_ _J_ L_ _L_ L_ _J_ _J_ n 227 

v^i3 Viz Viz Viz Viz Viz "~ 



z, L_ 1 1 1 1 1 1 n 99^ 

Viz Viz Viz Viz Viz Viz Viz u — 

a 1 111111 



0.223 



V13 V13 V13 V13 V13 V13 V13 

7 4- ±& 0.223 

Q _J_ _J_ _J_ _j_ _j_ _J_ 1_ 

v 7 !! Viz Viz Viz Viz Viz Viz 

9 0^ 0.219 



0.222 



V2 

10 -J- -i- -±- -±- -±- -±- 213 

JT ^ 0.195 
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Table 4: The total crystal-field splitting AE^ of p configuration (with Is coupling neglected) and 

(2) 

the average absolute values of E n in the crystal- field potentials given in Table 1. All the values are 
given in M2 units. 



— 

CF 



No. 



n 



AE^ 



1 

2 
3 
4 
5 



0.385 
0.381 
0.374 
0.369 
0.368 



ct 2 = v/2/25 = 0.283, A£^[ = 0.600, A£ ] 



0.600 
0.656 
0.683 
0.692 
0.693 



0.267 
0.262 
0.250 
0.239 
0.231 



horn 



max 



0.693 (Eqs 6-9) 
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Table 5: The total crystal-field splittings AE^ and AE^ of d 1 configuration (with Is coupling 



neglected) and the average absolute values of E. 



(2) 



and E ( n ] 



in the crystal-field potentials given in 



Tables 1 and 2. All the values are given in M2 and M4 units, respectively 



— 

CF 



No. 



H 



AE^ 



1 

2 
3 
4 
5 



0.385 
0.381 
0.374 
0.369 
0.368 



0.572 
0.572 
0.572 
0.572 
0.572 



0.229 
0.227 
0.221 
0.217 
0.213 



O 2 = V 2 / 35 = 


0.239, A£^ B 


= 0.488, ASlH = 0.676, AE^L = 0.756 (Eqs 6-9) 










av 


aeW 




a.v 



6 

7 

8 

9 

10 

11 

12 

13 

14 

15 



0.287 
0.280 
0.277 
0.276 
0.273 
0.269 
0.266 
0.265 
0.261 
0.251 



0.476 
0.363 
0.449 
0.437 
0.463 
0.426 
0.398 
0.549 
0.555 
0.564 



0.152 
0.174 
0.169 
0.169 
0.164 
0.159 
0.159 
0.137 
0.132 
0.113 



CT4 = 72/63 = 0.178, A£~Vf n = 0.363, A£^ m - 0.503, AS^L = 0.564 (Eqs 6-9) 
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Table 6: The total crystal-field splittings AE^ , AE^ and AE^ of f 1 configuration (with Is coupling 

neglected) and the average absolute values of E^ , E^f 1 and E„^ in the crystal-field potentials given 
in Tables 1, 2 and 3. All the values are given in M2, M4 and Mq units, respectively. 



No. 


H {2) 
n CF 


a.v 


ab( 2 ) 




av 


1 

2 
3 
4 
5 




0.385 
0.381 
0.374 
0.369 
0.368 


0.600 
0.603 
0.607 
0.608 
0.608 


0.190 
0.191 
0.192 
0.193 
0.193 




^4/75 = 


0.231, AS^L 


= 0.467, = 0.693, AS^L = 0.864 (Eqs 6-9) 










a.v 


AEW 


E (4) 
J-'n 


a.v 



6 
7 

8 
9 

10 
11 
12 
13 
14 
15 



0287 
0.280 
0.277 
0.276 
0.273 
0.269 
0.266 
0.265 
0.261 
0.251 



0.394 
0.417 
0.449 
0.458 
0.464 
0.478 
0.482 
0.399 
0.363 
0.358 



0.121 
0.119 
0.116 
0.117 
0.115 
0.113 
0.115 
0.129 
0.129 
0.131 



cr 4 = V 2 /" = 


0.142, A£ ( V B = 0.287, A^ m = 0.426, A£^ x = 0.531 (Eqs 6-9) 








n CF 


AE^ 

a.v 




a.v 



16 
17 
18 
19 
20 
21 
22 
23 
24 
25 
26 



0.239 
0.231 
0.228 
0.227 
0.225 
0.223 
0.223 
0.222 
0.219 
0.213 
0.195 



0.408 
0.326 
0.365 
0.379 
0.445 
0.420 
0.346 
0.422 
0.468 
0.481 
0.501 



0.107 
0.130 
0.120 
0.123 
0.106 
0.112 
0.127 
0.112 
0.097 
0.097 
0.072 



(J 6 = 



10 

13v^S3 



0.134, AS, 



0.271 



0.402 



0.501 (Eqs 6-9) 
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Table 7: The total crystal-field splitting AE& , AE^ and AE i& > of the 3 # 4 state and the average 

absolute values of En \ En^ and E^ in the crystal-field potentials given in Tables 1, 2 and 3. All the 
values are given in M2, M4 and Mq units, respectively. 



No. 

1 

2 
3 
4 
5 



n CF 

0.385 
0.381 
0.374 
0.369 
0.368 



AEW 

0.504 
0.543 
0.562 
0.560 
0.524 



tin 

ai 

0.163 

0.170 

0.164 

0.164 

0.155 



cr 2 = 0.184, A£ 


^, = 0.370, = 0.570, A^ x = 0.781 (Eqs 6-9) 










A£( 4 ) 

a.v 




a.v 



6 
7 
8 
9 

10 
11 
12 
13 
14 
15 



0287 
0.280 
0.277 
0.276 
0.273 
0.269 
0.266 
0.265 
0.261 
0.251 



0.215 
0.227 
0.249 
0.241 
0.232 
0.230 
0.232 
0.231 
0.229 
0.196 



0.079 
0.073 
0.070 
0.071 
0.073 
0.064 
0.067 
0.077 
0.076 
0.077 



horn 



cr 4 = 0.082, A£ ; 



0.165, AS, 



0.254. 



0.348 (Eqs 6-9) 





w (6) 

n CF 




A£( 6 ) 


E ® 

J-'n 


a.v 


16 


0.23 


f 


0.202 


0.05 




17 


0.231 


0.192 


0.069 


18 


0.228 


0.249 


0.051 


19 


0.227 


0.212 


0.060 


20 


0.225 


0.224 


0.062 


21 


0.223 


0.233 


0.059 


22 


0.223 


0.245 


0.053 


23 


0.222 


0.233 


0.061 


24 


0.219 


0.258 


0.052 


25 


0.213 


0.208 


0.059 


26 


0.195 


0.206 


0.058 



horn 



a 6 = 0.071, AS, 



0.143, AS, 



0.220, AS, 



0.301 (Eqs 6-9) 



18 



FIGURE CAPTIONS: 

Fig.l. Crystal-field splitting of |«7= 1) state - geometrical interpretation (x is the energy) a) general 
case: AS^ = x\- x 2 , b) A£^ n = x\- x 2 , and c) A£^ m = A£m2x = x 2 - x 3 . 

Fig. 2. Nominally allowed A£^ (bold solid borders) and the actual AE^ (dashed borders) ranges 
of the total splittings of the H4 state subjected to the iso-modular ~H CF . The A£^ om are also given 
(thin solid lines). 
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20 



AE * 

M, 



0.781 



0.7 - 



0.5- 



0.570 0.562 

0.504 



0.370' 



0.3- 



0.1- 



AE * 



0.4 -- 



0.3 



0.2-- 



0.1-- 



k = 2 



0.348 



0.254 



0.165' 



AE f 



0.4 -- 
0.3- 

0.249 

0.196 0.2 
0.1 + 



k = 4 



0.301 
0.220 

0.143 



0.258 
0.192 



k = 6 



Fig.2 
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